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We analyzed the pharmacokinetics of C3G on data from twelve subjects, after 2-week multiple dosing 
of black bean (Phaseolus vulgaris, Cheongjakong-3-ho) seed coat extract, using the mixed effect analysis 
method (NONMEM, Ver. 6.2), as well as the conventional non-compartmental method. We also 
examined the safety and tolerability. The PK analysis used plasma concentrations of the C3G on day 
1 and 14. There was no observed accumulation of C3G after 2-week multiple dosing of black bean 
seed coat extract. The typical point estimates of PK were CL (clearance)=3,420 1/h, V (volume)=7,280 
L, Ka (absorption constant)=9.94 h -1 , ALAG (lag time)=0.217 h. The black bean seed coat extract was 
well tolerated and there were no serious adverse events. In this study, we confirmed that a significant 
amount of C3G was absorbed in human after given the black bean seed coat extract. 
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INTRODUCTION 

Antioxidants are molecules capable of inhibiting the oxi- 
dation of other molecules. Oxidation is a chemical reaction 
that transfers electrons from a substance to an oxidizing 
agent. Oxidation reactions can produce free radicals and 
these radicals can start chain reactions. When the chain 
reaction occurs in a cell, it can cause cell damage or death. 
Antioxidants may terminate these chain reactions by re- 
moving free radical intermediates, and inhibit other oxida- 
tion reactions [1], 

Hence, plants and animals maintain complex systems of 
multiple types of antioxidants. Either low levels of anti- 
oxidants or inhibition of antioxidant enzymes may cause 
oxidative stress and may damage or kill cells. As oxidative 
stress appears to be an important initiator of many human 
diseases, the use of antioxidants has been intensively stud- 
ied, particularly as treatments for stroke and neuro- 
degenerative diseases [2,3]. Antioxidants have also been 
widely used as dietary supplements and have been inves- 
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tigated for their potential in the prevention of diseases, 
such as cancer and ischemic heart disease [4]. 

In the black bean seed coat, there are functional in- 
gredients, such as anthocyanins, which are known as pow- 
erful antioxidants (0.87 ~ 23.52 mg/g), as well as dietary fi- 
bers (80%). The black pigmentation is due to accumulation 
of anthocyanins in the epidermis palisade layer of the seed 
coat and cyanidin-3-glucoside (C3G, 80.9% of total content) 
is the most abundant anthocyanine [5]. 

There are some reports on the antioxidative effect of an- 
thocyanins extracted from the black soybean seed coat [2-4]. 
However, only limited information is available on the an- 
thocyanin's kinetics when the black beans were given to 
humans. In this study, we analyzed the pharmacokinetics 
(PK) of C3G from black bean using the mixed effect analysis 
method, as well as the conventional non-compartmental 
method. The safety and tolerability of black bean seed coat 
extract, when given 1 g per day orally for 2 weeks, have 
also been examined. 



METHODS 

Inclusion and exclusion criteria 

Healthy volunteers, aged 20 to 45 years, with no clin- 
ically relevant conditions identified based on medical his- 
tory, physical examination, laboratory tests and electro- 
cardiography (ECG), were eligible for inclusion. Subjects 
with any history that could result in alteration of C3G me- 
tabolism or with hypersensitivity to beans were excluded. 

ABBREVIATIONS: C3G, cyanidin-3-glucoside; PK, pharmacokinetics; 
ECG, electrocardiography. 
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Table 1. Subject demographics* 



Number of subjects 
Sex 

Age (years) 
Height (cm) 
Weight (kg) 



12 

Male: 6, Female: 6 

29.5 (24-44) 
166 (158-185) 

59.6 (47.0-78.4) 



"Continuous variables are shown as median (range). 



Finally, 12 subjects were enrolled (Table 1). 

Study design 

An open-label, randomized, multiple dosing clinical study 
was conducted at a single center (Clinical Trial Center, 
Seoul St. Mary's Hospital). All subjects were given 1 g of 
black bean seed coat extract every morning for 14 days. 
The trial was designed and monitored in accordance with 
the good clinical practice guidelines of Korea and with the 
principles of the Declaration of Helsinki. The independent 
institutional review board (Seoul St. Mary's Hospital) ap- 
proved the protocol before execution of the trial, and all 
participants gave written informed consent. 

Study procedures 

1. Subject control 

Subjects were admitted to the Seoul St. Mary's Hospital 
clinical trial center on the evening of day -1 and day 13. 
During the outpatient period, subjects were required to 
maintain low-phenolic diet and to report daily to the clinical 
trial center for dosing and diet. They were dosed at approx- 
imately 9 AM, prior to breakfast, and were instructed to 
fast for at least 10 hours prior to dosing on days 1 and 
14, and for at least 4 hours thereafter. 

2. Blood sampling 

On the morning of days 1 and 14, subjects had 7 ml of 
peripheral venous blood sampled predose and at 0.25, 0.5, 
1, 1.5, 2, 4 and 6 hours after dosing. Blood was collected 
in heparin tubes placed on ice, and then centrifuged at 
3000rpm for 10 minutes at 4°C. The separated plasma was 
stored at — 70°C until assayed for anthocyanins. 



AB SCIEX, Foster City, CA, USA). C3G and malvidin-3-gal- 
actoside (M3G), as the internal standard, were purchased 
from Extrasynthese (Genay, France). 

Plasma samples were acidified immediately with 6 M 
HC1 of 1/20 volume after centrifugation at 3,000 rpm for 
10 min at 4°C and frozen at — 70°C until analysis. For anal- 
ysis, plasma samples (500 lA) were prepared by diluting 
250 tA 0.1% formic acid, after spiking 25 lA internal stand- 
ard at 800 ng/ml. HLB cartridges (30 mg, 1 cc, Waters) were 
previously conditioned with 1 ml MeOH and 1 ml 0.1% for- 
mic acid; then, samples were loaded and cartridges were 
washed with 1 ml 0.1% formic acid. Finally, the analytes 
were eluted with 500 /A 80% ACN with 0.1% formic acid. 
The fractions were evaporated after collection and the resi- 
due was reconstituted in 200 iA of 10% MeOH with 0.1% 
formic acid, and 10 lA was injected into a LC/MS/MS sys- 
tem [6]. 

In the HPLC system, Luna C18 column (5 urn, 100x2.0 
mm, Phenomenex) was used with a gradient mobile phase 
consisting of 0.1% formic acid:methanol with 0.1% formic 
acid (80 : 20 to 20 : 80, v/v) and a flow rate of 500 uVmin. 
The MS/MS system was operated in the positive ionization 
mode with electrospray and multiple reaction monitoring 
mode. Transition ions at m/z 449.18 to 287.30 and 493.16 
to 331.04 were monitored for C3G and M3G, respectively. 

The lower limit of quantification for C3G was 0.2 ng/ml, 
with calibration curves ranging from 0.2 to 50 ng/ml. The 
precision (relative standard deviation) and mean intra- and 
inter-day accuracies were < 10.00% and 99.36 to 100.7%, 
respectively. These results demonstrated that the plasma 
concentration analysis was reliable over the given range. 

Non-compartmental analysis 

Initially, the pharmacokinetic analysis of C3G was per- 
formed using non-compartmental analysis (NCA) (Phoenix 
WinNonLin 6.2; Pharsight, Mountain View, CA, USA). The 
actual blood sampling times were used in the pharmacoki- 
netic analysis. The maximum concentration (Cmax) of drug 
in serum at Day 1 and Day 14, and the time to reach C ma x 
(T m ax) were obtained from the observed values. The area 
under the plasma-concentration-time curve until the last 
observation (AUCiast) was estimated using the log-linear 
trapezoidal rule. The terminal half life was calculated by 
dividing 0.693 by A. e . 



3. Safety 

Adverse events, physical examinations, vital signs, labo- 
ratory tests (clinical biochemistry, hematology and uri- 
nalysis) and 12-lead ECG were assessed during the study. 

Anthocyanin extraction procedure 

Twenty kg of black beans (Cheongjakong-3-ho) havested 
in 2010 was used to extract anthocyanin at a certified GMP 
unit. The bean seed coat was filtered and extracted for 24 
hours with 80% ethanol. The extract was finally packaged 
into 7,400 capsules containing 200 mg extract in each 
capsule. 

Assay of the plasma concentration of C3G 

Plasma concentration of C3G was determined by a vali- 
dated method using high-performance liquid chromatog- 
raphy coupled with tandem mass spectrometry (API4000, 



Statistical analysis 

Statistical analysis was carried out using R software 
Version 2.11.1 (R Foundation for Statistical Computing, 
Vienna, Austria). Descriptive statistics were used to de- 
scribe the PK parameters of C3G and the paired t-test was 
used to evaluate the changes in PK parameters on day 1 
and day 14. 

Mixed effect PK model development 

The concentration-time data for C3G were analyzed using 
NONMEM version 6.2 (Icon Development Solutions, Ellicott 
City, MD, USA) with the G77 FORTRAN compiler. The 
first-order conditional estimation (FOCE) method with in- 
teraction was used throughout the model building. 

One- and two-compartment open models with first-order 
elimination were tested to estimate the clearance (CL), cen- 
tral volume of distribution (VI), peripheral volume of dis- 
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tribution (V2), intercompartmental clearance (Q), absorp- 
tion constant (Ka) and lag time (ALAG), using the AD VAN 
subroutines. Models were selected based upon a decrease 
in the objective function value (OFV) of more than 3.84 
(p-value 0.05 in approximate X 2 distribution) and improve- 
ment in individual plots, as well as scatterplots. A log nor- 
mal distribution was assumed for inter-individual varia- 
bility ( v), and the PK parameters of the / h subject (Pj) were 
therefore described using the following equation: 



Model evaluation 

Ninety-five percent confidence intervals (CIs) for mean 
population PK parameters were determined using the 
re- sampling technique based on the bootstrap method. One 
thousand re-sampled datasets were estimated using the 
population PK model. The model was also evaluated by vis- 
ual predictive checks (VPC) with 1,000 simulated datasets. 



p. = TVP ■ exp( Vj ) 

Where TVP represents the typical population value of PK 
parameters, such as CL, V, Ka and ALAG. The inter-in- 
dividual variability eta ( v) for each PK parameter was as- 
sumed to follow a Gaussian distribution with a mean 0 and 
a variance co 2 . Possible correlations between the inter-in- 
dividual variabilities for CL and V were also estimated. The 
inter-occasion variabilities (IOV) for CL, V, and Ka, i.e. the 
variability within one individual between study days, was 
also examined with the exponential random effects term: 



RESULTS 

Non-compartmental analysis 

Paired t-tests were used to evaluate the changes in PK 
parameters on day 1 and day 14. The results are shown 
in Table 2. 

There were no significant changes in AUCi as t and C ma x 
on Day 1 and Day 14. No accumulation of C3G after 2-week 
multiple dosing of black bean seed coat extract was found. 



P jk = TVP-exp(n j + K jk ) 

Where Pjk is the f 1 individual parameter value at the day 
(occasion) of k that differs from the typical individual value 
by an additional random effect Kjk • Vj ■ Kjk and were as- 
sumed to be symmetrically distributed with a mean of 0 
and a variance of co 2 . 

As for the residual error, the additive, proportional and 
combined forms were tested. An example of the combined 
error form is shown, as follows: 

Y = IPRED- . + Js ~ . 2 + (e . .xIPMED..) 2 

Where IPREDy is the individual predicted concentration 
and Yy is the measured concentration of the / h individual 
at the * sampling time, and e,j is residual errors. Residual 
errors (e) include intra-individual variability, assay error 
and model misspecification. They were also assumed to fol- 
low a Gaussian distribution with a mean 0 and variance 
a\ 

Age, height, weight, sex and creatine clearance were 
screened as potential covariates of the parameters using 
Generalized Additive Modeling (GAM) implemented in 
Xpose version 4.0.4. 

In forward selection of covariates, variables that de- 
creased the minimized OFV greater than 3.84 (p<0.05) and 
decreased the omega value were selected. Covariates that 
did not increase the minimized OFV more than 3.84 (p< 
0.05) in backward elimination were removed from the 
model. 



Table 2. Pharmacokinetic parameters (Non-compartmental an- 
alysis) 



Parameter 


Day 1 


Day 14 


p-value* 


AUCi as , (ng • h/ml) 


4.26±1.40 


4.85±1.65 


0.0737 


C max (ng/ml) 


2.10±0.60 


2.17±0.68 


0.6659 


T„ al (h) 


0.73±0.29 


0.75±0.34 


0.8837 


Tl, 2 (h) 


1.58±0.33 


1.47±0.29 


0.3056 



*Paired t-test. 



Final model 

A one-compartment first-order elimination model with 
combined residual errors was chosen as the basic model. 
After establishment of the basic model, the inter-occasional 
variability between day 1 and day 14 was tested for all PK 
parameters, but was finally found relevant only for absorp- 
tion rate constant (Ka). Thereafter, the effects of covariates 
on relevant PK parameters were explored. There was no 
covariate that decreased the OFV and inter-individual 
variability. The final model was as follows. 

^ ) = BoseX C K (e" ** - e K * ) 

a e 

Where K a is absorption rate constant and Ke is terminal 
elimination constant. 

K aj = TVK a Xexp( Vj +K jk ) 



Table 3. Final model parameter estimates* 
Parameter Esitimate %RSE ^ 00 /,!!^ aP , 



CL/F 3,420 1/h 8.68 3,390 (2,890-4,070.3) 

V/F 7,280 1 6.84 7,250 (6,429.5-8,300) 

Ka 9.94 IT 1 30.7 9.6 (5.49-20.3) 

ALAG 0.217 h 3.67 0.216 (0.19-0.23) 
Between-subject 
variability 

(«CL 28.9% 42.0 26.5 (12.4-37.8) 

oj v 21.3% 49.2 20.0 (8.1-29.5) 

fflKa 102.5% 51.0 81.4 (41-130) 

pCL~V 0.883 0.945 (0.686-1.002) 
Residual error 

cr p „ p 21.3% 9.30 20.8 (17.8-24.6) 



*CL/F, oral clearance; V/F, volume of distribution; Ka, absorp- 
tion rate constant; ALAG, lag time; RSE, relative standard 
error; p, correlation coefficient; Prop, proportional. **95% CIs 
obtained from estimation of 1,000 bootstrap-resampled datasets. 
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Where K a j represents absorption rate constant of the / ' 
individual, TVK a is typical value of the absorption rate con- 
stant and Kjk is the interoccasional variability of f h 
individual. 

The population PK parameter estimates are summarized 
in Table 3. Basic goodness-of-fit plots for final PK model 
are presented in Fig. 1, 

Model evaluation 

Ninety-five percent confidence intervals (CIs) for mean 
population PK parameters were determined using the 
re-sampling technique based on the bootstrap method. One 
thousand re-sampled datasets were estimated using the fi- 
nal population PK model. The 2.5th and 97.5th percentiles 
of the mean population PK parameters were regarded as 
the lower and upper 95% CI limits, respectively (Table 1). 




0.0 0 5 1 0 1.5 2.0 0.0 0.5 1.0 1.5 2.0 2.5 3.0 

Population Prediction (ng/mL) Individual Prediction (ng/mL) 




0.0 0.5 1.0 1.5 2.0 2.5 3.0 0 1 2 3 4 5 6 

Individual Prediction (ng/mL) Time after Dose (h) 



Fig. 1. Goodness-of-fit plots for the final population PK model of 
C3G. Black line, line of identity; gray line, LOESS (locally weighted 
regression smooth line). IWRES, individual weighted residuals. 



VPCs were performed by using simulated concentrations 
of 1,000 virtual datasets (nsub=l,000 in the $SIMULATION 
block, 1,000 virtual patients) from the final model. The re- 
sult from VPC showed that the PK model gave acceptable 
predictive performance. Curves for the 5th, 50th and 95th 
percentiles of concentrations were overlaid on the observed 
concentrations on day 1 and day 14 (Fig. 2). 

Safety 

Overall, black bean seed coat extract was well tolerated. 
All 12 subjects completed the study. 

There were no serious adverse events (AEs) or discontin- 
uations due to AEs during the entire study period, but only 
one subject reported one minor adverse event, i.e. acute ton- 
silitis, which was mild and unlikely related to the con- 
sumption of black bean seed coat extract. The review of the 
ECG data did not reveal any safety concerns. There were 
no clinically relevant changes in terms of physical examina- 
tions, ECGs, vital signs or laboratory tests, including clin- 
ical biochemistry, hematology and urinalysis. 

DISCUSSION 

In this study, we evaluated the pharmacokinetics, as well 
as safety and tolerability, of C3G after 2-week admin- 
istration of black bean seed coat extract. To predict the ex- 
posure to C3G in healthy subjects population, simulation 
based upon the mixed effect analysis method was also used. 
There were no severe adverse events during the entire 
study period but only one minor adverse event, which was 
unlikely to be related to the consumption of black bean. 

Our data demonstrate that anthocyanins, especially C3G, 
are present in sufficient quantities in plasma for the PK 
parameters to be determined after oral intake of 1 g of the 
extract. The elimination of anthocyanins in plasma ap- 
peared to follow first-order kinetics, consistent with results 
from other studies [7,8]. 

In this study, the mean half-life of C3G (calculated with 
the CL and Vd estimated from the model) was 1.5 h, which 
was more than one hour shorter than that of other studies 
[8,9]. The differences in the formulations used in the stud- 
ies may have caused this difference. In our study, we used 
black bean seed coat extract, which was highly soluble and 



DAY 1 DAY 14 




Fig. 2. Visual predictive check plots 
of the final PK model; 1,000 virtual 
patients (1,000 datasets) were simu- 
lated using the final PK parameter 



Time after Doss (h) Time after Dose fh) estimates. 
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the absorption rate was very fast when given fasted; in con- 
trast, in other studies, subjects ate the raw fruit pulp, pow- 
der or other extract [8,9]. With its short half- life in human 
plasma, there was no observed accumulation of C3G after 
2-week multiple dosing in our study. 

In terms of the AUC of C3G, black bean seed coat extract 
lg was equivalent to about 250 ml of acai pulp and to about 
30 g of black raspberry powder [8,9]. 

As well as C3G itself, its metabolites produced in the 
body or bowel and other anthocyanins in the black bean 
seed coat may also have antioxidative effects. Therefore, 
further studies are needed to compare the antioxidative ef- 
fects of these natural products. 
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